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The possibility that fundamental discreteness implicit in a quantum gravity theory may 
act as a natural regulator for ultraviolet singularities arising in quantum field theory has 
been intensively studied. Here, along the same expectations, we investigate whether a non¬ 
standard representation, called polymer representation can smooth away the large amount 
of negative energy that afflicts the Hamiltonians of higher-order time derivative theories; 
rendering the theory unstable when interactions come into play. We focus on the fourth- 
order Pais-Uhlenbeck model which can be reexpressed as the sum of two decoupled harmonic 
oscillators one producing positive energy and the other negative energy. As expected, the 
Schrodinger quantization of such model leads to the stability problem or to negative norm 
states called ghosts. Within the framework of polymer quantization we show the existence 
of new regions where the Hamiltonian can be defined well bounded from below. 
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I. INTRODUCTION 

The standard model of particles has foundations on local quantum held theories having opera¬ 
tors of mass dimension d < 4. These operators are justified in order to implement the requirements 
of stability and unitarity without further elaborations [1]. However, when going to higher energies 
it is commonly believed that higher-order operators will play a key role in describing fundamental 
physics. This may be particularly true when they involve higher time derivatives since the new 
modes that arise allow to describe effects from a high scale. Usually these new modes are very high 

‘Electronic mail: creyes@ubiobio.cl 



2 


when the coupling of the higher-order time derivatives are suppressed by the high scale as taken 
in effective approaches. Higher-order operators have received increased attention over the years. 


They have been investigated in the context of loop quantum gravity 


tion 


0Q> 


corrections 


causality and stability [14], fine-tuning 


-Q, 


243] , Lorentz symmetry viola- 


21ri23] and nonminimal couplings 


154181], the hierarchy problem [HJ, 20], radiative 
The presence of higher-order derivatives in 


the gravitational sector are a key ingredient in order to achieve a consistent renormalization in the 


semiclassical approach, where matter fields are quantized over classical curved back grou nd 
Higher time derivative theories were introduced long time ago by Ostrogradsky 


H- 


311 ]. The ap¬ 


proach is based on a variational formalism and involves a Lagrangian L(x, x, x ..., x^ = Cjfr) 
and an extended Hamiltonian H(x, 7 t±,7Tx ■ ■ ■, 7r x („)) of 2 n variables. Ostrogradsky showed that the 
formalism leads a higher-order Hamiltonian not bounded from below as can be seen in the non¬ 
quadratic momenta terms that appear. This is the classical Ostrogradsky instability of higher-order 
time derivative theories which can be avoided in a few cases, for instance when the models has 
constraints 32j. 


The quantization of higher time derivative theories can be implemented introducing a change 
of variables in order to transform the original Lagrangian into a sum of decoupled normal-order 
Lagrangians. In general one of these Lagrangians has large negative energy leading to the instability 
or alternatively, by changing the vacuum state, to an indefinite metric theory [33|. The instability 
of the Hamiltonian has received much attention and has been tackled from different perspectives, 


such as phase space reduction 


Q, 


and Euclidean-path reduced amplitudes 


3443 . 6 ], complex_canonical transformations H, VT symmetry Q- 
4l|, gravitational ghosts and tachyons [42|, |43[. In 


quantum field theory, Lee and Wick showed that by imposing the negative norm states to decay it 
is possible to preserve unitarity order by order in perturbation theory J44, 4fj|. Resurgence of such 
ideas have been used to solve the hierarchy problem [19, 20] and in higher-order effective models 


with Lorentz symmetry violation 


46-48]]. 


In this work we study the stability of higher-order time derivative models within the framework 
of polymer quantization. In particular, we focus on the Pais-Uhlenbeck (P-U) model. The polymer 
representation is a non-standard representation of quantum mechanics inspired by some results 
that emerge from loop quantum gravity. The possibility of space discreteness that appears in loop 
quantum gravity Q has served to improve the convergence of quantum field theories jfj, 7] and 


cosmological singularities [501 ]. In this paper our main goal is to test whether the fundamental 
discreteness implicit in the polymer representation allows to improve the stability of higher-order 
theories. The polymer quantization has been considered in several studies such as two-point func- 
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tions 
namics 


58], compact stars 


5 


-53], cosmology 5C5^], central forces [55], 56[^ 


and low energy limits 


60, 


righer space derivatives [53], thermody- 


6l| 


The organization of this work is as follows. In section [IT] we introduce the P-U model and we 
explicitly show the origin of the instability in the Schrodinger quantization. In the third section we 
give a basic review of the polymer formalism. In section IIVI we polymer quantize the P-U model 
and solve the Hamiltonian eigenvalue equation. From the previous results we analyze stability in 
the region of validity of the theory. In the last section we give the conclusions. 

II. THE PAIS-UHLENBECK MODEL 


The P-U model consists of an harmonic oscillator coupled to a higher-order term described by 
the Lagrangian 


T 1 -2 1 7 2 9-2 

L = —mx - kx - x , 

2 2 2 

where g is a small parameter. The equation of motion is 

, d 2 x d A x 

_ fa _ m __ 9 _ = 0 . 

Inserting the plane wave ansatz x(t) = xoe~ iyt produces the four solutions 


(1) 


( 2 ) 


y = ±i 


I m ± \Jm 2 — 4 kg 


2 g 


(3) 


We define the two positive solutions according to 


Wi = U 0 ] 


l- 

2s 


<X>2 = W 0 1 


'l + x /TW4F 

2s 


(4) 


where we introduce uq = \Jk/m as the usual frequency and g = ^ with s a small dimensionless 
parameter. Now, in the limit e —>• 0 the first solution tends to the usual harmonic solution uj\ —>• ujq 

and the second one ui 2 —> blows up. One expects this behavior of the last solution in theories 

v £ 

with higher-order time derivative theories indicating a possible window to physics at higher energy 
scales. In order to avoid imaginary solutions we impose s < s c = 1/4 where s c is the critical value 
at which the two solutions collapse uj\ = W 2 - 

The conjugate momenta to x and x are defined by the expressions 

dL d fdL\ 

(5) 

( 6 ) 


^ dx dt \dx) ’ 


7r = 


dL 

dx 
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From ©, they are given by 


p = mx + gx 
7r = —gx 


(3) 


(7) 

( 8 ) 


The Hamiltonian is constructed via the extended Legendre transformation H = px + ttx — L which 
after substitution yields 


2 • 2 
_ 7 r mx 


kx 2 

2 9 2 +^ + — 


Let us consider the new set of variables 

lo\x -it/g 


Xl = 


2 2 ’ 


X2 = 


ujx - ir/g 

9 9 ’ 

^2 - Wf 


and 


Pi=p-(m- gul)x , p 2 =p-(m- gu\)x , 


(9) 


( 10 ) 


( 11 ) 


where m! = y m 2 — 4 kg. They allow to define the ladder variables 


oi = 


a 2 = 


m'uj i / 1 

xi + i——pi 

Zn \ mail 


m'uj2 / . 1 

X2 + l— - P2 


( 12 ) 


2 h \ m'uj2 

Using these variables the Hamiltonian can be expressed as the sum of two decoupled harmonic 


oscillators 


tvjj\ , + + . hui 9 , + + 

H = —{aia\ + a[ai) -— [a 2 a ' 2 + a^a 2 ) , 


(13) 


where we refer to the first and the second term as the positive and negative sectors of the theory. 

The quantization of such model follows by imposing the usual commutation relations in the 
extended phase space [x,p\ = ih and [x, k] = ih. With the use of the canonical variables [xi, p\ ] = 
ih and [x 2 , j5 2 ] = ih one can check that the creation and annihilation operators satisfy 


hi, (i ] 


= 1 , 


a,2, a\ 


= 1 . 


(14) 


To find the ground state wave function denoted by 4io, consider the explicit action of the operators 
in the Hilbert space 


xip = xip , 
Xl/j = xi/j , 


Pip = -9t— 
ox 

* / • f, 

TTW = 

OX 


(15) 

(16) 
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Inserting these expression into the previous creation and annihilation operators (I12|) we identify 
two realizations for the vacuum. The first one is to define the vacuum as the one annihilated by 
ai and a 2 


fil*0 = “2*0 = 0, 


(17) 


leading to 


= A/~ 0 e 9(u ™“ 2) ( ±2 +“i“2 x2 )+t uiu2xi 


(18) 


where A/"o is a normalization factor. The Hamiltonian turns out to be 

H = foui(JVi + -) - fuj 2 (N 2 + ^) , (19) 

where N\ = a\a\ and N 2 = a\a 2 are the number operators of positive and negative particles 
respectively. We see that the energy is not bounded from below since one can always create more 
negative particles. 

A different vacuum Tq amounts to change the annihilation operator in the negative sector and 
to maintain the previous in the positive sector, namely 


Ol*0 = 02*0 = 0 

with a' 0 = at. The vacuum state 'hf, in this case is 


Tg = A/q 


( 20 ) 


( 21 ) 


Note that Tq can be obtained in (1181) performing the transformation uj 2 —>• — uj 2 . The Hamiltonian 
is found to be 

H' = huJi(Ni + — ) + frw 2 (N' 2 + i (22) 

where now [a' 2 , ] = —1 and N 2 = —a 2 a 2 is the new number operator. In this case the theory has 

positive defined Hamiltonian but the price to pay is to end up with negative norm states or ghosts 
that may threaten the conservation of unitarity and with non-normalizable wave functions. 


III. POLYMER REPRESENTATION 

In quantum mechanics when dealing with the adjoint operators x and p usually one encounters 
some technical problems due to their unboundedness. Therefore, it is convenient to switch to the 
exponentiated versions U(a ) and V(/3 ) 

U(a) = e iai , V(P) = e^ h , (23) 
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whose action are defined by 


U(a)ijj(x) = e iax ^(x) , V(P)il>(x) = ip(x + /3) 


(24) 


for all state V’( x ’) in the Hilbert space L 2 (R). Both operators U(a) and V(/3) satisfy the Weyl- 
Heisenberg algebra 

U(a)U(at) = U(a + a') , 

vw)v(^) = vw + p '), 

U(a)V(/3) = e- iaP V(P)U(a) , (25) 

where the parameters a and /3 have (length ) -1 and length dimensions respectively. From the 
above algebra one can obtain the usual commutations relations [x,p] = ih. Due to the Stone-von- 
Neumann theorem any representation of the commutation relations have the form of the operators 
(1231) . modulo unitarity transformation, since the two operators U(a ) and V(/3) are strongly con¬ 
tinuous in their parameters, see [49j, l62j and references therein. 

In the polymeric construction one starts with a graph given by a countable set of points in the 


real line, denoted by 7 = {xj : j E N}, with some requirements 49]. We define the functions 
associated to a graph 7 as 


/(*) = 


fj x = x i 

0 x 0 7 

and their Fourier transform functions f(k) given by 


(26) 


—ixjk 

? £ 


(27) 


satisfying the relation 


Ei fi 


2 x? n < 00 for n = 0 , 1 , 2 ,... 


(28) 


We denote by Cyl 7 the space of all cylindrical functions f(k ) and Cyl the union of all Cyl 7 over 
all graphs 7 . We add to the space Cyl all the limits of Cauchy sequences, that is the Cauchy 
completion which is called the polymeric Hilbert space denoted by H po i y endowed with the scalar 
product 


(e 


—ikxj I —ikx 


j ) = 5 : 


I'Ej * 


(29) 


Recall, this is an alternative form to view the Hilbert space within the construction of the rigged 


Hilbert space D C H CO*, see Ref. 


63] 
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The main differences between the polymeric representation and the usual of quantum mechanics 
is that Hp 0 iy is a non-separable space and has an intrinsic fundamental discreteness leading to a 


nonequivalent representation, see Ref. [§4j]. To be precise, the action of the operators U(a) and 
V(/3) given in Eq. (|24D is well-defined in H po i y , however, in the polymer representation there is no 
self-adjoint operator p = such that the second equality in (1231) is satisfied, that is to say, the 

momentum operator p is not well defined on H po i y . This is due to the fact that V(/3) is not weakly 
continuous in the parameter f3, as can be verified using the modified product with the Kronecker 
delta in Eq. (12911 . Nevertheless, one can approximate the operator p with the expression 


„ ih 

p = - 

A* o 


- V(-p 0 


(30) 


where po is a fundamental length scale associated with a possible discreteness of space, coming 
from a more fundamental theory. The above approximation is natural, at least in the distributional 
sense, since if we take the limit as po —>• 0 we recover the usual momentum operator in L 2 (M). 

IV. STABILITY AND HIGHER-ORDER TIME DERIVATIVES 


In section m we have expressed the P-U Hamiltonian as two decoupled harmonic oscillators, for 
instance H = H\ — H 2 . Their quantum counterparts represent normal particles and nonstandard 


ones producing negative energy sometimes called Lee-Wick particles 
variables (fTOl) and (fTTl) we find 


44 


45|. Using the new set of 


pi 


Hl = \ hX " + 2 m' 

H 2 = \ k * xl + ^ 


(31) 


where kj = muj 2 - with j = 1,2. In other words, the P-U model involves two oscillators with 
the same mass m' and so taking advantage of this fact we polymer quantize the system as two 
individual harmonic oscillators. 

The polymer Hilbert space H po i y = H po i y (x 1 ) (g> H po i y (x 2 ) comprises the polymeric spaces for 
each oscillator. In the Hilbert space we have the action of the operators 


Z7i(ai)V>(aq) 
ViWi)^(xi) 
U2(a 2 )'ip(x2) 
V 2 (/ 3 2 )^ (x 2 ) 


_ e zaixi 


V>(xi) , 


= e 


^{x 1 + /3i) , 
ia2X ^(x 2 ) , 
V>( X 2 + P 2 ) ■ 


( 32 ) 
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Considering the wave function ij>(x 1 , 3 : 2 ) G H po i y we arrive at the Schrodinger equation 


^ + 2m' 


1 2 1 2 1 2 \ / 

Pi ~ ^ k 2 x 2 - ^P2 I nZl,X 2 ) 


= Eip(x i,x 2 ) 


(33) 


where E is the total energy of the system and the momentum operators are 

Pi = -^(Ci( W /2)-Vi(- W /2)) , 

P2 = (v2(m 2 /2) - F 2 (-/i 2 /2)) , (34) 

with and ^2 the fundamental lengths associated to H po i y (x 1 ) and H po i y (x 2 ). 

With the ansatz ■(/’(in,^ 2 ) = ^(xi)^^) and considering the cylindrical function for each 
oscillator j, namely 

= ^2tpj{x jie )e~^’ ek , (35) 

l 

together with Eqs (1341) . we obtain the equations for the coefficients 


+ Pj) + - Pj) 

/ 2E jP 2 i 


= 2 - 


3*0 + rt X lt 


frwjdj 






(36) 


where we have introduced E = E\ — E 2 and cfc = , . The previous equations suggest that one can 

J lib U/J 

hnd solutions supported at the uniformly spaced points Xj/ = Xo +i^j f° r some xo G [0, minj/n,}). 
Indeed, given the parameters of the equations {ip(xj/)}(2 = _ 00 , by using (1361) one can construct a 
solution i/jj supported on the lattice r yi 0}Pj := { Xj/ € M | x 3 j = xo + £/Xj, £ € Z}. Using this in Eq. 
(1361) yields 


^j(xj,e+ 1 ) + ipj(xj\/-i) 

^ f 2 E jP) ^(xo + tpj ) 2 

y frwjd* d',j 


4>j{xj,i) ■ 


(37) 


Without loss of generality, let us consider both graphs based on the point xo = 0. Hence, with 
Xj/ = ifij, ipj(kj) = YtL-oo ipj( x j,e)e~ lkjeflj and multiplying the equation (1371) by and 

summing over £, we arrive at 


1 - 


%i 

Hujjdj 




2cos(kj/j,j)'ifij(kj) = 2 


(38) 










9 


where kj € Normalizing and arranging the terms we get 

j fij fij 

( E ■ d 2 \ 

il>"{kj) + 2 d 2 j -i- + -^(cos (kjHj) - 1) ipj{kj) = 0 . (39) 

ynuij Hj J 

Finally, by making the change of variables Zj = — ^ +7r G (0, n) we write 

i/’jizj) + [aj ~ 2qj cos(2^)] ipj(zj) = 0 , (40) 

where * - 4Aj‘, a, = ^ (# - l) and X 3 = »/d, is a dimensionless parameter. Elation 
(1401) is the well-known Mathieu equation in its canonical form. We seek for periodic solutions of 
the Mathieu equations, since ipj( 0 ) = V ; j( 7r ) by construction. 

In order to approximate perturbatively to the harmonic oscillator we take Ai to be small, which 
produces the the following asymptotic expansion for the first oscillator 

a i,n.(Ai) = — 8 Aj 4 + 4(2n + 1)A^ 2 — — (2n 2 + 2n + 1) 

+0 (A 2 ) , (41) 


for n = 0, 1, 2,.... Considering the asymptotic expansion for a\ we can write 

Ei,n = -— ^(2ra + 1) — — (2n 2 + 2n + 1)A 2 ^ + 0{\\) . (42) 

For the second oscillator we have two alternatives: the first one is to consider /ro d 2 which 

analogously produces 

E' 2 ,m = - — ^(2 m + 1) — — (2m 2 + 2m + 1)A|^ + 0(A|) . (43) 


The second alternative is to consider a large A 2 leading to the expansion 

E 2 ,m = 


^ + hw 2 ^-m 2 + 0{ A 2 6 ) . 


(44) 


With exception of the ground state we have that the negative energy of the system increases 
without limit. This case may be interesting to analyze from a field point of view, since in this case 

51]. For the rigid rotator case in which Ai 


the propagators turn to be suppressed by the high scale 


is large as well as A 2 , achieved for example taking large values of the original mass parameter m, 
one can see that the energy spectrum goes as E n m ~ |(A 2 n 2 — A 2 m 2 ) which can be more negative 
with respect to the Schrodinger representation for the same occupation numbers (n, m). In this 
case there is no improvement for stability. Let us focus on the case Ai <S 1 and A 2 <S 1 for which 
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FIG. 1: Comparison of the P-U energies as a function of m evaluated at n in the Schrodinger and polymeric 
representation. 



FIG. 2: The polymeric P-U energy E n ^ m as a function of n and ra¬ 
the total energy E n ^ m of the system can be written as 

E n ,m = -— + 1) — — (2n 2 + 2n + 1)A 2 ^ 

- ^ ^(2m + l)-^(2 m 2 + 2m + l)AlV (45) 

Analogously to the Schrodinger quantization the high energy oscillator seems to lead to the in¬ 
stability, however we will show below that in certain regions the Hamiltonian can be defined well 
bounded from below. It is important to emphasize that in the absence of operators connecting 
the two Hilbert spaces the negative energy is not to serious and the problem of instability appears 
upon introducing the interactions. 

Let us consider the constraints imposed on the number of normal particles n that follows from 
the absence of any polymeric effect in quantum mechanics. From the first term in (14511 . it can be 
seen that the positive-oscillator corrections become significant or 0(1) at the value n = \fh\. For 
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example considering the vibrational modes of a carbon monoxide molecule with mass m = 10~ 26 
= 10 15 s _1 , and estimating the polymer scale to be = 10~ 19 m, we find 
In addition, in several approaches to polymer quantum mechanics a cutoff in the 
energy eigenvalues has been justified for the harmonic oscillator in order to implement a consistent 
renormalization [<3(]]. This upper limit is key to define our effective region, since now, in contrast 
to what happens in the usual P-U model for higher values of the occupation number m and for 
fixed n = n the polymeric energy is bounded from below, see Fig. [I] In the general setting where 
both occupation numbers n and m vary freely the total energy of the system is well bounded from 
below, avoiding possible stability problems, as shown in Fig. [2j 


kg and frequency 
n « 10 7 , see 1491. 


V. CONCLUSIONS 

In this work we have analyzed the stability of theories containing higher-order time derivatives 
within the framework of polymeric quantization. For this we have focused on the well-known 
Pais-Uhlenbeck model with fourth-order time derivatives in the Lagrangian. Using a canonical 
transformation we have cast the theory into a sum of two decoupled harmonic oscillators, one with 
positive energy and the other with negative energy. The negative-energy oscillator is responsible 
for the instabilities that arise in the presence of interactions. 

We have shown that the discrete nature of the polymer Hilbert space introduces corrections 
in the energy spectrum of the P-U model which allows to define a region of positive defined 
Hamiltonian. For this we have set a cutoff from observational constraints from quantum mechanics 
due to the absence of any polymeric effect and further motivated by a consistent renormalization 
program. We have established an effective region defined by small values of the parameters Ai 
and A 2 at which the theory has a well bounded Hamiltonian. However, for the case with large Ai 
and A 2 , we have found that the instability shows up for very low occupation numbers with no 
improvement with respect to the usual Schrodinger quantization. We leave for future investigations 
the inclusion of interactions and the case of large A 2 in the context of quantum field theory. 
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